One Dimensional Magnetized TG Gas Properties in an External Magnetic Field 
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With Girardeau's Fermi-Bose mapping, we have constructed the eigenstates of a TG gas in an 
external magnetic field. When the number of bosons N is commensurate with the number of 
potential cycles M, the probability of this TG gas in the ground state is bigger than the TG gas 
raised by Girardeau in 1960. Through the comparison of properties between this TG gas and 
Fermi gas, we find that the following issues are always of the same: their average value of particle's 
coordinate and potential energy, system's total momentum, single-particle density and the pair 
distribution function. But the reduced single-particle matrices and their momentum distributions 
between them are different. 
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I. INTRODUCTION 

The problem of a one-dimensional (ID) hard-core gas 
was first studied classically by Tonks [l| . Girardeau 0, @| 
continued the research in quantum, and put forward the 
Fermi-Bose mapping to solve energy spectrum and the 
corresponding wave function of a one dimensional Tonks- 
Girardeau (TG) gas. After that, Lieb and Liniger Q 
considered N Bose particles interacting via a repulsive 
5- function potential with a coupling constant 7. In fact, 
TG gas can be considered as the Lieb-Liniger gas when 
7 — > 00. 

During the past decade, the study on TG gas has un- 
dergone a rapid development in theoretical and experi- 
mental p5l-[T7l]. In particular, A. Lenard [j| indicated that 
the off-diagonal parts of the one-body density matrix and 
the momentum distribution show distinct differences in 
both the TG gas and the Fermi gas. G. J. Lapeyre et 
al. [§[ studied the momentum distribution of a harmon- 
ically trapped gas. The most important progress is the 
realization of the strongly interacting bosons [l6| in a 
one-dimensional optical-lattice trap and the TG gas [jjj 
by trapping 87 Rb atoms by trapping them with a combi- 
nation of two light traps. 

We have known some properties of one dimensional 
TG gas, and we are interested in its properties and fea- 
tures in external magnetic field. Consequently, we first 
reveal some properties of both gases in the same exter- 
nal magnetic field and compared them. With the Fermi- 
Bose mapping @, 0|, we obtained all the eigenenergies 
and eigenstates of this TG gas. When the number of 
bosons N is commensurate with the number of potential 
cycles M, the TG system in an external magnetic field 
will be in the ground state with a bigger probability than 
the TG gas raised by Girardeau in 1960. Many properties 
of the TG gas and the Fermi gas are always of the same 
even if it's time-dependent, such as their average value of 
particle's coordinate and potential energy, system's total 



momentum, single-particle density and the pair distri- 
bution function. But both their reduced single-particle 
matrices and momentum distributions are different. 



II. TG GAS IN AN EXTERNAL MAGNETIC 
FIELD 

In a TG gas, the boson is assumed to have an "impen- 
etrable" hard core characterized by a radius of a. From 
Girardeau's work, with the hard core radius a — > 0, the 
interparticle interaction is given by 
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Such an interparticle interaction could be represented 
by the following subsidiary condition on the wave func- 
tion ip : 

tl>(x ir -- ,x N ,t)=Q if x t = x j ,l^i<j^ N. (2) 

In order to let the TG gas wave function, which must 
be symmetric with respect to the permutations of any Xi 
and Xj, satisfy this condition ©, we need the ideas of 
Fermi-Bose mapping proposed by Girardeau [1, 0] . That 
is, using the fact that the Fermi wave function, denoted 
by ip F , satisfies condition ^ naturally and is antisym- 
metric, a bosonic wave function vb B of a TG gas can be 
constructed by 



in which 



,x N ,t) = A(xi, ■ ■ ■ ,x N )^ F (xi,- ■ ■ ,x N ,t) 



N 



A(x lr ■ ■ ,x N ) = ^Qsgn(a 



(3) 
(4) 



and 
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As the ground state of a bose system is non-negative [H 
mapping (|3|) for the stationary ground state of a system 
reduces to a simplified form 



ip B {xi, ■ ■ ■ ,x N ) = \ip F {xi,- ■ ■ ,x N )\ 



(5) 



To illustrate our general ideas, we now study the case 
that magnetized bosons in an external magnetic field 
B(x) = —Bcos(2ljx). For simplicity, we take B{x) inde- 
pendent of t. therefore, 



H = 



N 
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i=l 



h 2 d 2 

2m dx 2 



in which 



V(x) 



fj,B cos(2oja;) 
0, 



V(xi 



< Xi < £, 
else. 



(6) 



(7) 



We suppose L = Mir, where M is an integer and fi is 
the magnetic moment of atom. For the sake of simplicity, 
we assume that both N and M are odd (As a result of 
the introduction of A(x 1 , • • • , x N ), tb^ is periodic if N is 
odd and antiperiodic if N is even [2. H3|)- From Bloch 
theory, the corresponding energy spectrum of system in 
periodic potentials has the energy band structure. We 
use the notation n as the band index and k as the Bloch 
wave vector in the first Brillouin zone. For convenience, 
we use a = {n, k} to denote both the band index and 
the Bloch wave vector. Not considering the interparticle 
interaction, every particle in x G [0, L/u] is governed by 



h 2 d 2 

2m dx 2 



Substitute z = 



\±B cos(2cja;) 



ip a (x) = E a cp a (x). (8) 



UJX, q 



m[iB 2mE , 

—z — j and A = — ^ — ^ into equa- 



tion (|5|), then this equation becomes 

d 2 <p(z) 



dz 2 



+ [A - 2qcos(2z)]ip(z) = 0, 



z G [0,L]. 



(9) 



It is called Mathieu's Differential Equation, which can 
be solved with Randall B. Shirts's [18|] method. Using 
Bloch's theorem, we first set 



<p(z) = exp(ivz)u(z) 
ip(z) = ip(z + n), 
u(z) — u(z + L), 
21 

v = — 
M 



I G 0,±1,--- ,± 



M - 1 



(10) 

(11) 
(12) 

(13) 



Since u(z) is periodic with period 7r, it can be expanded 
in Fourier series: 



<p(z) — exp(ivz) c n exp(i2nz). 



(14) 



Then substitute 1)140 into ©, we obtain an infinite 
symmetric tridiagonal matrix equation, that is, eigen 
equation. 
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By truncating the matrix in each direction (centered at 
the smallest diagonal element v 2 ) at sufficiently large di- 
mensions, approximations to the desired eigenvalues can 
be obtained to any desired precision. We could obtain 
the eigenvalue and eigenfunction easily using Mathemat- 
ica or Matlab. In practice, the eigenvalue and eigenfunc- 
tion converge quickly as the dimension increases when 
q is not too large. For example, if we choose v = 6/7 
and q = 1, the relative difference of the lowest eigenvalue 
between truncating the matrix at 21-D and 2001-D is 
only 10 -38 . For large values of q and for high orders, we 
can use asymptotic expansions instead (2ll . |23 | . since the 
matrix dimension that needed to get accurate eigenvalues 
and eigenfunctions becomes large. 

With Girardeau's Fermi-Bose mapping ((3]), the eigen- 
states of the hard-core Bose system is given by the Slater 
determinant 



•0 B (a; 1 , • • • ,x N ) = 



-^(^1 ) ' * ' ) % N ) 



(fa^X,) (fa^X,) ■ 
<f ai (x 2 ) fa 2 (x 2 ) ■ 

f ai (x N ) <Pa 2 (x N ) ■ 

The total energy of the system is 

N 



Va N (#i) 

V? ajv (x 2 ) 
V a „{x N ) 



(16) 



(17) 



As a result of the Bose-Fermi mapping, the energy spec- 
trum of the Bose and corresponding Fermi system are 
identical. When the temperature is 0, the system will 
be on the ground state — N particles on the N lowest 
eigenstates, respectively. 
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Calculations show that the energy gap AE between the 
hrst and the second Bloch bands decreases with the de- 
crease of M, and when M is large, the energy gap AE is 
independent of M. The energy gap between the first and 
the second Bloch bands as a function of magnetic field B 
and frequency u are plotted in Fig. 1, respectively. At 
5 = 0, raised by Girardeau 0, the energy gap of TG 
gas is the smallest (not 0). The energy gap increases as 
B and u increase. As^the Boltzman Distribution Law 

notes, — — oc exp( — 7777)1 with AE = Ek — Ey, n k and 

n , are particle numbers on the corresponding states. So 
the system have a bigger probability in the ground state 
with large B and u> than the case B = when the to- 
tal number N — M. That is to say, one could realize 
ground-stated TG system more simple with large q and 
w than q = when the total number N = M. 



For the Fermi gas, the reduced single-particle density 
matrix can be simplified to 



AE (1(T 23 J) 4 




o)(10'm _1 ) 



FIG. 1: Plots of the energy gap between the first and the 

second Bloch bands as a function of magnetic induction B 

and frequency u. M = 9, m = 1.44 x 10~ 25 kg, p = 9.274 x 
10 _ 24j . T -i_ 



III. PROPERTIES OF THE TG GAS 
COMPARED WITH THE FERMI GAS 

In this section, we further study properties of the ID 
N magnetized hard-core Bosons in the external magnetic 
field 0, and then compared with the properties of ID 
N magnetized Fermions in the "same" state which hard- 
core Bosons occupy. 



A. Reduced single-particle density matrix and 
Single-particle density 

The reduced single-particle density matrix with nor- 
malization J p{x, x) dx — N is given by 

p B (x,x',t) = N J tp B *(x,x 2 , ■■■ ,x N ,t) 

x ip B (x' , x 2 , ■ ■ ■ ,x N ,t) dx 2 ■ ■ ■ dx N . (18) 



p F {x,x',t) =N J ip F (x,x a ,--- ,X N ,t) 

x ip F (x',x 2 ,- ■ ■ ,x N ,t)dx 2 --- dx N ( 19 ) 
= l P a (x,t)*f a (x\t). 



a— all states 
occupyed 
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FIG. 2: Density-plots of p(z, z') of the hard-core TG gas 
(a,b,c) and Fermi gas (a',b',c'), M = 7. Abscissa is z lt Ordi- 
nate is z 2 . (a,a) N = 3; (6,6') N = 5; (c,c') N = 7. 

Both the reduced single-particle matrices of hard-core 
bosons and fermions in the ground state in the exter- 
nal magnetized field ([7]) are plotted separately in Fig. 2. 
The multidimensional integral in equation (|18[) is evalu- 
ated numerically by Monte Carlo integration using Math- 
ematica. Because the relationship between x and z is 
proportional, for the sake of simplicity, we take z as a 
variable. From the figure we can see the off-diagonal ele- 
ments of the matrix of both the two kinds of gases decay 
quickly as N increases; The bright diagonal stripe of TG 
gas is thinner than Fermi gas under the same condition. 
This reflects the condensate properties of Bose gas; And 
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the shock of the off-diagonal elements p(z, M ir — z) of TG 
gas is smaller than Fermi gas. To see this more clearly, 
we have plotted the off-diagonal elements p(z, Mir — z) 
in Fig. 3. Sudden rise in the border is the result of using 
periodic boundary conditions. 
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FIG. 3: Plots of p(z,Mn - z), M = 7. p%(z,Mn - z) 
stands for p(z, Mir — z) of n hard-core Boson particles, so 
as p F (z, Mtt — z) 



Using the mapping ([3]), we can prove that the single- 
particle density p(x,t), normalized to N, of both hard- 
core Bosons and Fermions are equal, no matter whether 
they are on the ground-state or not. 



dx t 
dx„ 



(20) 



p B (x,t) = N J \ip B (x,x 2 ,--- ,x N ,t)\ 2 dx 2 
— N J \i) F (x,x 2 ,--- ,x N ,t)\ 2 Ax 2 

a— all states 
occupyed 

= p F (x,t) 



p(x,t) is just the element of the reduced single-particle 
matrix p(x,x',t) when x = x' . Fig. 2 shows that p(x) is 
cyclical, which indicates that particles tend to stay cycle 
at where its potential energy is low. 



B. Pair distribution function 

The pair distribution function, which is the probability 
of finding a second particle as a function of distance from 



an initial particle, normalized to N(N — I), is defined as 

D(x 1 ,x 2 ,t) 

= N(N-1) [ \ip B ( Xl ,x 2 ,--- ,x N ,t)\ 2 dx 3 --- dx N 



N(N-l) / \i> F ( Xl ,x 2 ,--- ,x N ,t)\ 2 dx 3 --- dx r 



9 ^2 \ ( P a (x 1 ,t)(p a ,(x 2 ,t) - ip a (x 2 ,t)ip a ,(x 1 ,t)\ 2 . 

(21) 



a, a =all states 
occupied 



From the above equation we can see that the pair distri- 
bution functions of both hard-core Boson gas and Fermi 
gas are identical. In fact, the average of particle's coor- 
dinate and potential of both the two gases are identical, 
too. Because they involve absolute values of the wave 
functions only. Also we can say this feature is deter- 
mined by the fact that the single particle density of both 
hard-core Boson gas and Fermi gas are exactly of the 
same. For the ground state of a system with N particles 
governed by Hamilton (|6|), the pair distribution function 
is 



D(x 1 , x 2 



N-l 

E 

a.ct'—0 



\<Pc.( X l)<P a > (X) - PafaJV* / (^l)l 2 



Figure 4 shows density-plots of the pair distribution 
function of the ground state for different number of par- 
ticles, when x 1 = x 2 , D{x 1 ,x 2 ) — 0. It is the result 
of the impenetrable hard-core intcrparticlc interaction. 
Just like the single particle density, the pair distribution 
function reflects the same periodicity and particles tend- 
ing to tarry at the low potential energy. As the number of 
particle increases, the black strip becomes thinner. One 
of the reason is that the averaged distance between par- 
ticles reduces as the number of particles increases while 
the scope of the potential field is certain. In two distant 
regions, away from the diagonal in Fig. 4, D{z 1 ,z 2 ) is 
cyclical. This reflects the hard-core interaction between 
particles is local. 



C. Momentum distribution 

The momentum distribution, related to the reduced 
density matrix, is defined as 

n{k) = 1- [ p(x,x')e- ik ^- x,) dxdx', (23) 

27T J 

with the normalization 



n(k) dk = N. 



Actually, the momentum distribution is just the Fourier 
transformation of the reduced density matrix (|18[) . We 
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FIG. 4: Density-plots of the pair distribution function 
D(z 1 ,z 2 ). Abscissa is z x , Ordinate is z 2 . M = 7. (a) N = 2; 
(b) AT = 3; (c) JV = 5; (d) N = 7. 

know that the Fourier transform of the original function 
and the transformed function is one-to-one, as the differ- 
ence of the reduced density matrix p(x, x') of the hard- 
core Bose gas and Fcrmion gas, the momentum distri- 
bution is surely not the same. Yuan Lin and Biao Wu 
[ToT] have plotted the momentum distributions for both 
the TG gas and the free Fermi gas in a periodic Kronig- 
Penney potential, which shows the difference — Even 
a free Fermi gas has a broader momentum distribution 
than the most strongly interacting boson gas. We would 
like to add that, although their momentum distribution 
functions are different, their total momentum are of the 
same. With the equation 

(^ r + ^)^(«.-*J = 0, (24) 

we can prove 

N „ 

^2—A(x 1 ,---,x N )=0. (25) 



And then, go ahead, we can obtain 



P(t) B = P(t) F , where P = -ihj^ if-- (26) 



IV. SUMMARY AND CONCLUSIONS 

On solving the eigenfunction of the TG gas in the ex- 
ternal magnetic field, we found that the TG system with 
large B and u) has a bigger probability in the ground 
state than the TG gas raised by Girardeau in 1960 when 
the number of bosons N is commensurate with the num- 
ber of potential cycles M. It reveals that we can realize 
ground-stated TG system more easier. And then we have 
studied properties of magnetized TG gas and Fermi gas 
in an external magnetic field, comparing with each other 
when they are in the "same" state reveals that it is im- 
possible to distinguish them just from their average value 
of particle's coordinate and potential energy, system's 
otal momentum, single-particle density and pair distri- 
bution function. But we could distinguish them from 
their reduced single-particle matrices or their momentum 
distributions. In order to illustrate our results, we have 
plotted their reduced single-particle matrices and the mo- 
mentum distributions when they are in the ground state. 
Although their momentum distribution functions are dif- 
ferent, their total momentum are of the same. These 
results will be deeply helpful for understanding the TG 
gas. 
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